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ABSTRACT 

We investigate the relationship between the star formation rate (SFR) and dense molecular gas mass in the 
nuclei of galaxies. To do this, we utilize the observed 850 /im luminosity as a proxy for the infrared luminosity 
(Lir) and SFR, and correlate this with the observed CO (J=3-2) luminosity. We find tentative evidence that 
the LiR-CO (J=3-2) index is similar to the Kennicutt-Schmidt (KS) index (A^ « L5) in the central ^L7 kpc 
of galaxies, and flattens to a roughly linear index when including emission from the entire galaxy. This result 
may imply that the volumetric Schmidt relation is the underlying driver behind the observed SFR-dense gas 
correlations, and provides tentative confirmation for recent numerical models. While the data exclude the 
possibility of a constant Lir-CO (J=3-2) index for both galaxy nuclei and global measurements at the ~80% 
confidence level, the considerable error bars cannot preclude alternative interpretations. 
Subject headings: stars: formation - galaxies: ISM, evolution, stai^burst - radio lines: ISM, galaxies 



1. INTRODUCTION 

The rate at which stars form in galaxies has long been pa- 
rameterized in t erms of a power-law involving the gas density 
(ISchmidllll959h . HI and CO observations of galaxies have 
shown that the star formation rate (S FR) is related to t he sur- 
face gas density via Ssfr oc Y}^^^-^ (lKennicuttlll998h . The- 
oretically, a power-law index of ^L5 controlling the SFR is 
appealing. If a constant fraction of gas is converted into stars 
over a free fall time, the relation SFR oc pi.,^ results. 

Observations of molecular gas in galaxies have provided 
a more complex view of SFR relations. Local surveys have 
shown that the SFR (as traced by the Ljr) is proportional to 
the CO (J=l-0) luminosity (L'; Gao & Solomon 2004a) to the 
^ L5 p ower, consistent w ith observed surface density SFR re- 
lations (iKennicuttll 19981) . However, recent observations have 
revealed a tight, linear relation between HCN (J=l-0) and the 
Lir in galaxies dGao & Solomonir2004alf5 ). Because the J=l- 
rotational transition of HCN has a relatively high critical 
density (wcrit ^^-^10^ cm"-*), this has been interpreted as a more 
fundamental SFR relation such that the volumetric SFR is lin- 
early related to the dense, star-forming molecular gas. A lin- 
ear relation between Ljr and the luminosity from a high crit- 
ical density tracer CO (J=3-2) in a similar sample of galaxies 
has provided confirming evid ence for this relation a nd advo- 
cated a similar interpretation ( Narayanan et al.ll2005 l). 

An alternative interpretation to the Unear relation between 
Lir and dense gas tracers ha s been put forth from the the- 
oretical side. First, Krumhol z & Thompson! ( l2007h utilized 
models of individual giant molecular clouds (GMCs) with 
a lognormal density distribution function coupled with es- 
cape probability radiative transfer calculations. These authors 
found that the central issue driving the observed molecular 
SFR relations was the relationship between the fraction of the 
cloud's gas above the critical density of the molecular tran- 
sition. When, on average, the gas density is higher than the 
line's critical density (<ncioud> » «crit), the bulk of the gas is 
thermalized in the line, and the line luminosity increases lin- 
early with increasing mean gas density (<n>). In this scenario, 
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a relationship between SFR and Lmoi is expected to have index 
similar to the underlying (volumetric) Schmidt index. This 
is reminiscent of the observed relation between Lir and CO 
(J=l-0) in galaxies. Alternatively, when only a small fraction 
of the gas is thermalized {<n^\oud> « Wcrit), the line luminosity 
will increase superlinearly with <«>, and the SFR-Lmoi rela- 
tion will have index less than that of the underlying Schmidt 
in dex (e.g. the observed L ir-HCN J=1-0 relation). 

'Narayan an et all (l2007h furthered these models of GMCs 
by applying 3D non-LTE radiative transfer calculations to hy- 
drodynamic simulations of isolated galaxies and equal mass 
binary galaxy mergers. A key finding in these numerical mod- 
els was that the observed relations are only found when global 
measurements are made. Higher spatial resolution observa- 
tions of e.g. the nuclei of galaxies would probe gas in which 
a larger fraction of the gas is thermalized than in unresolved 
observations of the entire galaxy. In this case, <n> could po- 
tentially become comparable to the critical density of the line, 
and an SFR-Lmoi index of ^L5 would be expected even for 
high critical density molecular lines. In both sets of models, 
the fundamental relation is the Schmidt relation with index 
L5. The observed SFR-Lmoi relations were simply manifesta- 
tions of the underlying Schmidt law. 

A key difference exists be t ween the numerical mod- 
els of iKrumholz & ThompsonI (l2007h an d [Naravan an et alj 
('2007 1), and the interpre tations of Gao & Sol omonI (I2004 a'b'). 
.Narayanan et all (l2005h and IWu et al. (2005,). The models 
find the driving relation is the volumetric Schmidt relation 
whereas the latter set of observations cite a more funda- 
mental relation as that between SFR and dense molecular 
gas. Tests that distinguish between these interpretations are 
crucial to understandi ng global SFR relations in galaxies. 
[Narayanan et al.l (l2007h offered a direct prediction from their 
models that observations of molecular lines with critical den- 
sity higher than that of HCN (J=l-0) should break the linear 
trend seen between SFR and HCN (J=l-0). Specifically, for 
extremely high critical density lines, «critwill be even larger 
than <ngaiaxy> than in the case of HCN (J=l-0); In these cases, 
the relation between SFR and Lmoi is directly predicted to 
be sublinear. Indeed, confirming evidence for this trend has 
been found by Bussmann et al. (2008; submitted), who found 
a sublinear Lir-HCN ( J=3-2) index in remark able agreement 
with the predictions of lNaravanan et alj ( l2007h . 
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An alternative g e neric fea ture of both the 



D< 50Mpc 



iKrumholz & Thompson! (l2007h and iNarayanan eTaP (l2007h 
simulations which may serve as a test of the models is a 
superlinear SFR-Lmoi relation for high critical density tracers 
when <ngaiaxy> is sufficiently high. One potential mani- 
festation of this is a break in the Hnear SFR-HCN (J=l-0) 
relation in systems with extremely high infrared luminosity 
(e.g. hyper-luminous infrared galaxies; Lir> lO'^L©). Here, 
the models pre dict a steepening in the SFR -Lmni index 
toward A^=1.5 (Krumholz & ThompsonI l2007h. Tentative 
evidence for this may have been found by iGao et al.l (|2007t) 
in high redshift systems, though the potential contribution of 
active galactic nuclei (AGN) to the Ljr may drive a similar 
signature. 

In order to avoid the muddying effects of central AGN, a 
potential alternative to this test is to observe the SFR-Lmoi re- 
lation in galactic nuclei, where the local mean gas density may 
be higher than the globally averaged mean gas density. In this 
Letter, we utilize literature data in order to examine the rela- 
tionship between SFR (traced by the 850/im flux) and dense 
molecular gas (traced by CO (J=3-2)). The aim is to help dis- 
tinguish between the competing interpretations of the linear 
SFR-dense gas relations in galaxies. 

2. LITERATURE DATA 

In order to investigate the relation between SFR and dense 
molecular gas over a variety of physical spatial extents in 
galaxies, we require a beam matched set of observations in 
both high critical density molecular line and SFR tracer with 
a sufficient number of galaxies. While matching the highest 
resoluti on and sensitivity infra red data of nearby galaxies (e.g 
SINGS; iKennicutt & etl2003b with high spatial re solution ob- 
serva tions of a high critical density tracer (e.g. iKrips et al] 
l2007h would be ideal, few galaxies exist at the intersection of 
such surveys. Most other traditional infrared surveys are in- 
sufficient as they typically report global measurements, rather 
than higher spatial resolution data. Further complications 
arise on the dense gas tracer side. Because of the prodigious 
observing time necessary to map large numbers of galaxies 
with an interferometer, most dense gas surveys are done with 
a single dish and report on lines at ^mm wavelengths (e.g. 
iBaan et al.ll2008l) . In order to obtain the highest spatial res- 
olution possible from single dish surveys, higher-frequency 
observations (e.g. A<1 mm) must be employed. 

The largest beam-matched sample of a dense gas tracer and 
LiR ti-acer is the JCMT CO (J=3-2) survey of Yao et al. (200l 
and th e SCUBA Local U niverse Galaxy Survey ( Dunne et al.l 
I2OOOI) . lYao et al.1 (l2003h used the 850/^m images of these 
galaxies to scale the total Lir to match the '^IS" beam of the 
CO (J=3-2) observations. As such, while the number statistics 
are relatively small, to our knowledge this sample comprises 
the largest available for this type of analysis. Moreover, using 
CO (J=3-2) as a tracer of dense molecular g as has the attrac- 
tive q uality that both global measurements jNaravanan et alj 
l2005h and simulated unresolved observations from simula- 
tions (Naravanan et al. 2007) have shown a linear relation be- 
tween Lir and CO (J=3-2) luminosity in galaxies. 

In order to avoid single galaxies at the extrema of the Lir 
range artificially biasing the fits, we impose nominal lumi- 
nosity cuts in the sample, considering galaxies in the range 
LiR=[10^, 2x 10" Lq]. This excludes galaxies which may be 
interacting, thus allowing us to use galaxy distance as a proxy 
for mean density (§[3]l. This results in a total sample size of 40 
galaxies, excluding only 4 galaxies from the parent sample. 
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Fig. 1 . — The LiR-CO (J=3-2) relation for galaxies in each of our distance 
bins. The slope is seen to flatten as observations probe more global measure- 
ments. 
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Fig. 2. — The LiR-CO (J=3-2) index as a function of limiting galaxy 
distance. The physical scale covered by the beam is on the top axis. Ob- 
servations of the nearby galaxies probe the nuclear region of the galaxies. 
In this high <n> regime, the Lir-CO (J=3-2) index tentatively lies near the 
volumetric Schmidt index, als o consis tent with recent mo deling efforts by 
IKrumholz & ThompsonI t2007l) and Naravan an et aljr2007l) . As observations 
become unresolved and more low density gas is folded into the beam, the in- 
dex tends toward unity. 



3. RESULTS 

We have binned our sample of galaxies into bins in distance. 
Because we are considering observations of a fixed beamsize, 
this is equivalent to binning the galaxies in terms of physi- 
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cal radial extent from the nucleus probed by the observations. 
The number of galaxies rises sharply to ^25 when consider- 
ing galaxies with a maximum distance of ~50 Mpc. We thus 
utilize 50 Mpc as the first distance bin in order to maximize 
the statistical significance of the result. 

In Figure [T] we show the Lir-CO (J=3-2) relation for each 
distance bin. In Fi gure [H we plot the be st fit log(LiR) - 
log(CO (J=3-2)) (U JOao & Solomonll2004al) slope as a func- 
tion of limiting distance. On the top axis we label the physical 
extent of the beamsize.The fits were done utilizing the publi- 
ca lly available Monte Carlo Markov Chain fitting packages 
of iKellyl (l2007h . The Bayesian routines assume the intrinsic 
distribution can be well approximated via a mixture of Gaus- 
sians, and returns the posterior probability distribution func- 
tion (pdf) of potential slopes to the fit. Unlike standard linear 
regression techniques, the fit is not done with the assumption 
that the abscissa values are known exactly while measurement 
errors exist only for the ordinate. Rather, errors are allowed 
in both axes. This is important a s errors are only reported for 
the literature CO measurements ("Y ao et afl2003h . and to con- 
form with standard literature fits, we are required to fit Lir as 
a function of the CO luminosity. The points in Figure |2] are 
the median of the returned distribution of slopes, and the error 
bars denote the standard deviation in the pdf. 

The slopes in the Lir-CO (J=3-2) relation tentatively show 
a similar trend to what would b e exp ected from the models o f 
iKrumholz & ThompsonI (l2007l) and iNarayanan et all (l2007h . 
Specifically, in the regions of high <n> where the bulk of 
the gas may be thermalized, the SFR-Lmoi relation for high 
critical density tracers tends toward 1.5 (e.g. the points at 
D«50 Mpc in Figure |2l this corresponds to the central ^1.7 
kpc for most of the galaxies in this bin). When considering 
more unresolved, global observations of galaxies which fold 
in significant amounts of diffuse gas (thus lowering the ef- 
fective <n> probed by the observations), the mean density of 
the galaxy drops below the critical density, and the observed 
SFR-Lmoi index is less than that of the underlying volumet- 
ric Schmidt index. Here, we see that the more unresolved 
observations approach a slope of unity, consistent with the 
linear slope found between Lir and CO (J=3-2) luminosity 
in a n early identical sample of galaxies by Narayanan et aT| 
(l2005h . We note that the slope is still moderately superlin- 
ear in these bins, th ough consistent with th e range of re- 
sults expected by the INarayanan et al.l (l2007h models. This 
occurs because the maximum distance bins still contain the 
relatively nearby (D<50 Mpc) galaxies, which increases the 
Lir-CO (J=3-2) index. Ideally one would like to exclude the 
nearby galaxies from the most distant bins, but low number 
statistics prevent this experiment. 

4. DISCUSSION 

The trends seen in Figure|2]may be consistent with a picture 
in which the observed SFR-Lmoi relations are driven globally 
by the relationship between <ngas> and ncrit- In contrast, if the 
Unear relation between Lir and HCN (J=l-0) or CO (J=3-2) 
were indicative of a more fundamental SFR relation in galax- 
ies in terms of dense molecular gas, one would expect the 
relation to remain linear even in the nuclei of galaxies. With 
a potentially model distinguishing result such as this one, ob- 
vious questions regarding its robustness arise. 

First, we caution that the error bars presented in the low- 
est distance bins in Figure |2] are rather large. This owes to 
the small sample sizes of these bins ('--^25 galaxies for the 
lowest distance bins). It is possible to investigate what the 



probability is that the same Lir-CO (J=3-2) index properly 
characterizes both the lowest and highest distance bins. To 
do this, we remind the reader that each point (and associated 
error bars) in Figure |2] is a pdf for potential Lir-CO (J=3-2) 
indices at each distance bin. In Figure [3] we plot the (nor- 
malized) difference in the pdf's from the highest and lowest 
distance bins in Figure |2] We additionally show the cumu- 
lative distribution function in the same plot. A value of in 
the pdf difference is expected at points when the same index 
characterizes both the high and low distance bins. As can 
be seen, the nuclear Lir-CO (J=3-2) index is systematically 
weighted toward larger numbers than the global value. That 
said, the error bars in the nuclear distance bin are not insignif- 
icant. The probability that the Lir-CO (J=3-2) indices from 
the nuclear and global distance bins are the same ±0.25(0.5) 
is ^-^17(33)%. Thus, while the tentative trends seen in Fig- 
ure |2] are indeed probable, they are by no means robust. In 
this sense, surveys to increase the number statistics will be 
required to confirm/refute this potential result. 

Second, we can inquire as to the validity of Lir as an SFR 
indicator. In particular, the Lir from galaxies has the poten- 
tial to be contaminated by central AGN. That said, the galax- 
ies in the sample presented here are relatively low l uminosity 
(global Lir < 10 L©). Both o bservational studies (|Tran & eti 
2001) and theoretical works (Chakrabart i et al suggest 
that galaxies with the Lir range investigated here have their 
IR luminosity largely powered by star formation. 

Third, we can question whether or not the trend in the Lir- 
CO (J=3-2) index seen with distance is equivalent to a trend 
in <n>. To do this, in FigureH] we plot the normalized <n> of 
the galaxy versus the galaxy distance. For the density deter- 
mination, the volume is derived by assuming the gas resides 
in a disk with radius equivalent to the physical extent probed 
by the beam at that distance, and a constant scale height in 
each galaxy. The molecular gas mass is taken from the CO 
(J=l-0) luminosity, utilizing a CO-H2 conversion factor ap- 
propriate for this sample of galaxies (lYao et al.l 120031) . The 
CO (J=l-0) observations are beam matched to the CO (J=3-2) 
observations presented here, so the molecular gas mass is in- 
deed the gas mass within the radial extent of the CO (J=3-2) 
beam. The mean density is normalized to account for a variety 
of uncertain parameters (e.g. disk scale height, molecular gas 
volume filling factor, density profile of GMCs). Because we 
are only concerned with the trend in <«gaiaxy> with distance, 
normalizing the density has no consequence. Figure |4] shows 
that the mean density probed by the molecular line observa- 
tions drops as we consider galaxies successively farther away. 
Alternatively said, unresolved observations of galaxies probe 
lower <«> than observations toward galactic nuclei. This im- 
plies that the tentative trend in Lir-CO (J=3-2) index seen in 
Figure|2]owes to the relationship between the mean density of 
the gas in the beam and nciit of the emission line. 

If larger samples verify the results of Figure |2] then these 
observed trends may hav e sever al im plications. The models of 
IKrumholz & ThompsonI (l2007l) and INarayanan eTaP (l2007l) 
which predict an SFR-Lmoi index of 1.5 when <n> is high 
are predicated on an underlying relation controlling the SFR: 
SFR oc «'^. The results presented here tentatively support a 
scenario in which the SFR can be described by such a power- 
law (rather than a linear relation), though we reiterate that the 
error bars in this study are substantial. 

In principle, we could perform the same experiment shown 
in Figure|2]with the CO (J=l-0) transition. Indeed, Nobeyama 
45m observations of CO (J=l-0) in many of the same galax- 
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Fig. 3.— Solid Line: Normalized difference in PDFs of Lir-CO (J=3-2) 
indices for lowest distance bin and lai'gest distance bin. Dotted line: Cu- 
mulative distribution function for shown PDF. This can be interpreted as the 
PDF of the LiR-CO (J=3-2) indices for the nuclear and global observations as 
being the same. 



a similar experiment with CO (J=l-0). Naively, one might 
imagine the Lir-CO (J=1-0) index to be '--^1.5 for all distance 
bins owing to the relatively low critical density of this ground 
state transition. However, the more likely scenario are results 
similar to the trend seen in Figure |2]- a slope of ~1.5 for 
the lowest distance bins, and a flattening to a slope of ^^1. 
The reason for this is global observations of galaxies in this 
relatively low Lir range have a linear Lir-CO (J=1-0) relation 
(iGao & Solomon 2004a b). It is only when considering galax- 
ies with LiR>10" Lq in the sample that unresolved observa- 
tions res ult in an Lir-CO (J=1-0) i ndex o f 1 . 5. In the interpre- 
tation of Ikrumholz & ThompsonI (l2007h and lNarayanan et all 
(l2007h . this would result from lower luminosity galaxies hav- 
ing most of their gas subthermal in this transition. 

Finally, we emphasize that though the results presented 
here advocate an interp retation of the linear L[R-dense gas 
relations consistent wifl i iKrumholz & ThompsonI (|2007|) and 
iNaravanan et alj (l2007h . this is not to say that some high crit- 
ical density tracers such as HCN (J=l-0) are not applicable as 
SFR tracers. By virtue of their linearity with Lir, insomuch 
that the Lir can be used as an adequate proxy for the SFR, so 
can global measurements of HCN (J=l-0) and CO (J=3-2). 
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Fig. 4. — The normalized <n> within the CO (J=3-2) beam as a function 
of galaxy distance. The <n> within the beam drops as a function of distance, 
showing that the more unresolved observations are probing lower <«> gas. 
This implies that the relationship seen in Figure[2]owes to a dropping <n> in 
the farthest objects. 

ies exist, thus providing th e opportunity fo r a beam-matched 
comparison to this dataset (lYao et al.ll2003h . However, the er- 
ror bars on these observations are sufficiently large that the 
slopes on the Lir-CO (J=1-0) indices for the lowest distance 
bins (D < 100 Mpc) are consistent with numerous interpreta- 
tions (including no correlation at all). 

The lack of clean data in low «crit transitions does not pre- 
clude us from conjecturing what results one might expect in 



5. CONCLUSIONS AND WAYS FORWARD 

We have studied the relationship between the SFR and 
dense molecular content in the nuclei of galaxies by utilizing 
the observed Lir as a proxy for the SFR, and CO (J=3-2) as 
a proxy for dense molecular gas. We find tentative evidence 
that the Lir-CO (J=3-2) index is superlinear in galactic nu- 
clei with value ^ 1.5, while flattening as the observations 
measure more global conditions. 

The error bars in this potential trend are considerable how- 
ever, and that the results here serve as motivation for future 
observational experiments investigating the LiR-dense gas re- 
lation in a variety of environments. One potential way to in- 
crease the sample size would be to obtain additional beam- 
matched sub-mm and CO (J=3-2) observations of the nuclei 
of nearby galaxies with 15 m sub-mm telescopes (e.g. the 
JCMT or ASTE). Alternatively, high spatial resolution obser- 
vations of nearby galaxies with ancillary IR maps (e.g. the 
SINGS sample), would provide a rich dataset with Lir and 
dense gas measurements for numerous regions with varying 
physical conditions in individual galaxies. Additional data to 
constrain the gas density (e.g. CO J=l-0) would allow the 
Lir-CO (J=3-2) index to be studied truly as a function of gas 
density, as opposed to galaxy distance. These experiments 
will help to confirm/refute the potential results of Figure |2] 
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